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WITH POROUS WALLS 

By William J. Nelson and Paul L . Klevatt 


SUMMARY 


A method of generating variable Mach number supersonic flow in a 
channel of fixed geometry by the removal of air through uniform porous 
walls is discussed. Calculated porosity distributions are presented 
for several minimum-length nozzles designed to operate at Mach numbers 
up to 2.0. The axial pressure gradient has been calculated for several 
constant-porosity walls over a range of Mach numbers. The applicability 
of these calculations to a two-dimensional tunnel is illustrated by 
comparison of calculated and experimentally determined pressure gradients 

in the 2-i- by Li -inch channel at Mach numbers in the range covered in the 
4 2 

present tests, 0.99 to 1.17. Schlieren photographs of the flow in the 
experimental channel are also presented. 


INTRODUCTION 


The development of variable Mach number supersonic wind tunnels 
has received a great deal of attention during recent years. Several 
schemes for varying the Mach number by changing the tunnel geometry 
have been tried. These include moving nozzle blocks (references 1, 2, 
and 3) and flexible walls (references 4 and 5 )• Another type of vari- 
able Mach number wind tunnel being investigated by the National Advisory 
Committee for Aeronautics and others utilizes a series of longitudinal 
slots in fixed tunnel walls. In these fixed- geometry tunnels the Mach 
number is varied by changing the pressure in the slots relative to the 
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the stream stagnation pressure. Preliminary results of this vork are 
available in references 6 to 9* 


The present investigation is a continuation of the slotted- tunnel 
vork but utilizes porous metallic vails rather than a finite number of 
passages as employed in the slotted tunnels. Either the slotted or the 
porous-vall tunnels offer the possibility of avoidance of transonic 
choking in addition to variation of the supersonic Mach number vith fixed 
geometry. Eesearch on both possibilities is in progress^ but the present 
paper is concerned only vith the problem of variable supersonic Mach 
number. This report contains a simplified analysis of the tunnel- vail 
characteristics required to generate shock-free supersonic flov in mini- 
mum distance together vith calculated and experimental pressure distri- 
butions along a constant-cross -section channel vith vails of knovn porosity. 
All calculations vere based upon the methods vhich are discussed in refer- 
ence 10 using tabular data from reference 11. The experimental data pre- 
sented vere obtained in a symmetrical channel of rectangular cross sec- 


tion ^ by 4 inches^. This vork vas done under steady-flov conditions 

vith air supplied at elevated pressures vhile atmospheric pressure vas 
maintained in the area surrounding the tunnel. The vork vas initiated 
at the Langley internal aerodynamics section in the summer of 19^8 and 
is continuing. 


SYMBOLS 


a sonic velocity 

g acceleration due to gravity 

p local static pressure 

Ap local pressure difference across porous surface 

M local Mach number 

Y velocity 

5 angle of flov relative to vail 

7 ratio of specific heats 

p mass density 

V cumulative sum of expansive deviations of stream from 

direction of flov at M = 1.0 

x distance along porous section 
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h 


height of channel 


Subscripts 


D 


design condition 


o 


stagnation condition 


n 


component normal to porous vail 


ANALYSIS 


The difference in Mach number between any two points in a two- 
dimensional supersonic tunnel bears a unique relationship to the cumu- 
lative sum of the angular deviations to which the stream has been 
subjected between those points. A complete discussion of the relation- 
ship between the Mach number and stream deviation can be found in 
reference 10, chapter 2. If the direction of flow at M = 1.00 is taken 
as a basis of reference, this relationship may be expressed as 


The generation of supersonic Mach numbers in conventional tunnels is 
accomplished by accelerating the air stream in a Laval type nozzle; the 
rate of acceleration and therefore the length of the tunnel being 
determined by the curvature of the walls. In a symmetrical, minimum- 
length, shock-free nozzle, one-half the design stream deviation Vp 

corresponding to the final Mach number Mp is produced by increasing 

the angle between the wall and the center line of the stream at a point; 
the remaining deviation, a result of the initial stream curvature along 
the opposite wall, produces a realinement of the flow to the tunnel 
center line. For the purposes of this analysis only the flow in the 
lower half of a symmetrical nozzle will be considered, thus, for the Mach 
number range under discussion, all characteristics with positive slope 
are defined as part of the first family expansion while those character- 
istics with negative slope are from the second family expansion, figure 1. 

Shock- free nozzles .- By initiating the first family expansion at a 
point and curving the wall downstream of this point to correspond to the 



4 


NACA RM L50B01 


direction of flow through, the second family expansion, the design Mach 
number will he reached in the minimum distance with continuous acceler- 
ation of the stream. The resulting nozzle is of the general type shown 
in the characteristic diagram of figure 1. Along the tunnel wall the 
velocity component normal to the tunnel axis increases from 0 to a maxi- 
mum at a value corresponding to v-j-,/2 at the point of initial turning; 

continued acceleration of the stream beyond this point is accomplished 
by curvature of the streamlines in the opposite direction with resulting 
reduction in the normal component to zero. The magnitude of these 
velocity components is shown for several minimum- length nozzles in 
figure 2(a). The curves, drawn as a function of v, are for completeness 
presented for both first and second family expansions although the region 
of practical interest is limited to that section of the curves drawn 
double weight. If the solid wall is replaced by a porous wall parallel 
to the plane of symmetry, the velocity component normal to the tunnel 
axis becomes the component normal to this wall; if the porosity of the 
wall is such that this flow may pass through the wall, minimum- length 
nozzle flow will be realized. The pressure variation as a function of v 
is also shown in figure 2(a). 

The rate of flow through any porous material is determined by the 
pressure difference across it and by its loss characteristics; thus, any 
flow distribution can be obtained by a suitable matching of the pressure 
differential across the material and its loss characteristics. The 
following analysis is based upon maintaining a constant-pressure chamber 
outside the porous wall and varying its resistance to flow to obtain the 
desired normal flow rate. Since the final pressure reached in the channel 
will be equal to that in the chamber outside the porous wall, the pressure 
increment across the wall at any point will be equal to the difference 
between the local static pressure and that corresponding to the design 
Mach number. This difference, for the shock-free nozzles of figure 2(a), 
is plotted in figure 2(b) . 

The velocity component normal to the wall and the pressure difference 


as a function of v in figure 3(a). As in the preceding figures, the 
first and second family expansions are differentiated by weight of line. 
For minimum-length nozzles, in which only the variation across the second 
family expansion is pertinent, the total change in the ordinate is rela- 
tively small until the higher Mach numbers are reached. To show the 
range of porosities that will be required for materials across which the 
pressure drop varies linearly with pg\T and corresponding curves 


across it have been combined into 



plotted 
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of 


p o a o 

Ap/Po 


and 



are also presented, figures 3 ("b ) and 3(c). 


As in figure 3(a), these curves are, for completeness, drawn for the 
first family expansion as well as for the second. From a comparison 
of these curves it is apparent that, for a given Mach number, a material 
for which the relationship between Ap and pgV' remains essentially 
constant offers a much closer approximation to the ideal porosity than 
one in which Ap varies with pgV^, since in the latter material the 
proportionality constant decreases to 0; this condition corresponds to 
a solid wall at the nozzle exit- 


Curves of the ratio of the weight flow component normal to the 
porous wall to the pressure drop across it, pgV n /Ap, as obtained from 
ficure 3(b), are presented in figure b as a function of distance along 
the nozzle for minimum- length nozzles designed to operate at Mach numbers 
of 1.1, 1.2, 1.5, and 2-0. For these curves T 0 was assumed constant 

at 500° F . From a consideration of the stagnation constants in the 
ordinates of figures 5(a) to 5(c), it can be seen that the porosity 
varies as the square root of the stagnation temperature. T e engt o 
nozzle required to generate a Mach number of 1.10 is approximately equal 
to one-half the tunnel height; increasing distances will be required to 
attain higher Mach numbers (x = h for M = 1*5 81111 x 2 for 

M = 2.0) . 


Corresponding curves based upon pressure drop as a function of 
pgV n 2 have not been prepared since it is apparent from figure 3(c) that 

shock-free flows will require a wide range of porosity. Although it is 
possible to attain good flow with such characteristics, it appears 
probable that over a range of Mach numbers materials across which the 
pressure drop varies linearly with the mass flow will provide the more 
satisfactory solution. 


Pressure distribution wi th constant porosity .- In the preceding 
section it was observed that the attainment of shock-free supersonic 
flow in a porous-walled channel requires a prescribed rate of air flow 
across the channel boundary. It was also shown that, for a minimum- 
length nozzle, the required local rate of air removal varied approxi- 
mately with the first power of the difference between the local pressure 
and the pressure corresponding to the final Mach number. To detemine 
the significance of departures from the calculated ideal porosity , the 
pressure gradient along the center line of a two-dimensional channel 
vith walls of constant porosity pgV n /Ap has been calculated. The 
significance of departures from the calculated ideal porosity is shown 
in figures 5(a) to 5(c) where the pressure gradient along the center 
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line of a two- dimensional channel has been determined for walls across 
which the weight flow to pressure difference at T Q = 500° was 0.0055* 

0.011, and 0.018. The pressures along a channel with walls of known 
porosity were determined by iteration. The normal component of the 
weight flow rate, pgV n , was determined for assumed values of V hy 

placing pgV n = p^Ia sin 8 and again from the known porosity hy multi- 
plying pgV n ^p hy the pressure difference corresponding to the assumed 
values of V. At equilibrium, the value of pgV n is independent of the 

method of calculation; the corresponding flow angle relative to the wall 
was then determined from the intersection of the two curves of pgV n 

against V. By repeating this process, the pressure distribution through- 
out the length of the channel was determined. 


The curves presented in figure 5 show that for — — £ = 0.0055 the 

Ap 

distance required to reach a Mach number of 1.10 is approximately one- 
half the tunnel height greater than that required for the minimum-length 
nozzle of figure b; for all higher Mach numbers the added distance 
required to reach equilibrium flow is much greater. Increasing the 
porosity to 0.011 (fig* 5(t>)) appreciably reduces the distance required 
to reach equilibrium for all Mach numbers; M = 1.10 is reached in 
approximately the same distance required by the minimum-length shock-free 
nozzle; greater than minimum length is required, however, for all greater 
Mach numbers and since the length of the tunnel would be determined by 
the maximum Mach number to be reached, this porosity also is too small. 


At a value of = 0 .018, figure 5(c), the initial expansion is 

Ap 

slightly excessive for Mach numbers below 1.20 but the extent of the 
over- expansion is small. For higher Mach numbers the distance required 
to reach equilibrium is somewhat in excess of the minimum although Mach 

x pgVji 

numbers up to 1 .5 are reached in ^ < 2 .0 . At higher values of 

the initial expansion at low Mach numbers is greater than that required 
for minimum length and subsequent compression is stronger; higher values 
of M can be attained, however, in very short nozzles. 


APPARATUS 


The general arrangement of the experimental apparatus showing a view 
of the test channel, the porous upper and lower walls, and the pressure 
tube used in making the static-pressure surveys is presented as figure 6. 
The contraction ratio of the entrance bell ahead of the test section was 
approximately 70:1. The 2^- by l^-inch rectangular cross section of the 
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channel was constant throughout the porous region. Detailed pictures 
of the porous walls are presented as figure 7. These surfaces were 
selected commercial grade sintered "bronze sheets consisting of ther- 
mally bonded bronze particles approximately 0.005 inch in diameter 
rolled to a thickness of 0.08 inch. The pieces used in this investigation 
had the smoothest surface finish of any samples readily available. The 
surface corresponded roughly to the surface of 120 wet-or-dry paper. 
Structural stiffness was provided by three longitudinal bars tapered to 

-!g-inch width at the point of contact with the sintered bronze; lateral 

support was provided at the ends only. 

]_ ' 

The inch -diameter pitot- static tube used in making the static- 

O 

pressure surveys had four static orifices 0 .031 inch in diameter; the 
total-pressure opening was 0.050 inch in diameter. Stagnation pressure 
and temperature were recorded at a point upstream of the entrance bell . 

A schlieren optical system, with 6- inch-diameter mirrors and a spark 
light source which produced intermittent flashes of approximately 
10 microseconds, was used in photographing the flow. 


EXPERIMENTAL EESULTS AND DISCUSSION 


The foregoing analysis has shown the desirability of a porous 
material that approximates a linear variation of pgV n with Ap for 
use in generating uniform supersonic flows. Such materials, however, 
were not available for use in this investigation. An experimental check 
of this method of calculation was therefore attempted using the material 
at hand and calculating the pressure distribution from an experimentally 
determined calibration. The air flow through the sintered bronze is 
shown as a function of the pressure difference across the material in 
figure 8. The apparatus used in the calibration is shown schematically 
in the upper part of this figure. In applying this calibration it has 
been assumed that the loss characteristics are a function of the velocity 
component normal to the surface and are completely independent of any 
component parallel to the surface. The calculated pressure gradients are 
compared in figure 9 with the experimental pressure distributions at 
supersonic Mach numbers up to 1.17; this limit was established by the 
pressure ratios available with existing equipment and does not represent 
a tunnel limitation. Corresponding schlieren photographs of the flow are 
presented in figure 10. 

In general, the experimental data are consistent with the calculated 
curves. Differences in absolute values of the pressure ratio as deter- 
mined from calculation and experiment are, at the forward end of the 
nozzle, attributed to failure of the calculations to allow for 
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non-uniformity of flow at this point. The abrupt turn assumed in the 
calculations can he realized only in the absence of boundary layer at 
the entrance section since the presence of a layer of subsonic air at 
the point of turning will appreciably reduce the rate at which the flow 
direction may be changed. Evidence of nonuniform flow at the tunnel 
inlet can be found in the schlieren photographs at the left side of 

figure 10, where multiple shocks ahead of the porous floor ^ = -0.4 

reach into the flow but do not completely span the channel . Experimen- 
tally determined velocity profiles at the inlet indicated that the 
displacement thickness of the boundary layer was approximately 0.02 inch 
for these tests . From one -dimensional considerations it can be shown 
that removal of this boundary layer alone will produce a Mach number of 
approximately 1.10 in the tunnel; the experimental curves would not, 
therefore, be expected to follow closely the calculated curves at the 
Mach numbers below 1.10. The reduction of the differences at the forward 
end of the channel between the experimental and calculated curves at 
M >1.1 is attributed both to the reduction in boundary-layer thickness 
and the increase in magnitude of the initial turn. The measured pressure 
ratios indicate that the flow at the center line of the main stream was 
supersonic ahead of the porous surface for the tests at M = 1.08 
and 1 . 17 ; this is ascribed to a thinning of the boundary layer ahead of 
the porous material . 

Maximum deviations in Mach number of d0 .01 are indicated by the 
experimental pressure distributions. Surface roughness is especially 
critical in this type of tunnel since supersonic velocities occur very 
close to the surface as a result of removal of substantial part of the 
air adjacent to the walls. The presence of the many Mach waves in the 
schlieren photographs is a result of the roughness of the sintered bronze 
floor. Parallelism of these waves gives some indication of the uniform 
velocity existing throughout the test regions. Normal shocks in the 
downstream section of the tunnel were caused by flow restriction at the 
end of the porous region where the walls were again solid. 


CONCLUSIONS 


The following conclusions are based upon the results of this pre- 
liminary investigation of the porosity requirements for a variable Mach 
number supersonic tunnel of fixed geometry: 

From analysis : 

(1) Shock-free supersonic flows can be obtained in a tube of 
constant cross section by removal of air through porous walls to a 
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constant -pressure chamber "by proper distribution of porosity along 
the tunnel . 

(2) The porosity requirements for shock-free nozzles vary with 
distance along the channel, Mach number, and with stagnation tem- 
perature . 

(3) Flows possessing only small compressions can be obtained 
in a constant-cross-section channel by removal of air across walls 
of constant porosity to a chamber of constant pressure. 

From experiment: 

(1) The measured pressure distributions are consistent with 
the calculated results. 

(2) The rate of expansion at the tunnel inlet is appreciably 
reduced by the presence of the boundary layer ahead of the porous 
material . 

(3) Flows with Mach number variations of iO.Ol along the center 
line of a fixed-geometry tunnel were obtained throughout the speed 
range of the tests (M = 0.99 to 1.17) by removal of air through 
commercial grade sintered bronze walls. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Xangley Air Force Base, Ya. 
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(a) v n/ a o and P/Po against V. 

Figure 2.- Flow conditions for shock-free nozzles. 
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(b) Ap/po against V. 


Figure 2.- Concluded. 
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(a) — against V. 
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Figure 3 »- Porosity requirement for shock-free flow. 
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Figure 3-- Continued. 
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Figure 3«- Concluded. 
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Figure 4.- Required porosity distribution for minimum-length shock-free nozzles. T Q = 500° R. 
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Figure 5.- Theoretical center-line pressure distributions for several Mach numbers 

constant pV n g/Ap. T 0 = 500° R. 
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Figure 5*- Continued. 
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Figure 5*- Concluded. 
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Figure 6.- General arrangement of experimental setup 
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(a) Top view of siritered-bronze floor. 



(b) Photomicrogram of sintered floor enlarged 25 diameters. 



(c) Bottom view of sintered-bronze floor showing lateral supports. 


Figure J.- Porous sintered-bronze floor. 
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Figure 8.- Porosity calibration of sintered bronze used in walls of experimental channel. 
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Distance along porous section, x/h 




Figure 9.- Pressure distribution along center line of 2^ - by ^ -inch 

porous-wall tunnel. 
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Figure 10.- Schlieren photographs of flow along 2 ^- by 4^ -inch porous- 

Inac/T 


wall tunnel. 
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